Diagnostics for the ground state phase of a spin-2 Bose-Einstein condensate 
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We propose a method to determine the singlet-pair energy of a spin-2 Bose-Einstein condensate 
(BEC). By preparing the initial populations in the magnetic sublevels and ±2 with appropriate 
relative phases, we can obtain the coefficient of the spin singlet-pair term from the spin exchange 
dynamics. This method is suitable for hyperfine states with short lifetimes, since only the initial 
change in the population of each magnetic sublevel is needed. This method therefore enables the 
determination of the ground state phase of a spin-2 87 Rb BEC at zero magnetic field, which is 
considered to lie in the immediate vicinity of the boundary between the antiferromagnetic and 
cyclic phases. We also show that the initial state in which relative phases are controlled can be 
prepared by Raman processes. 

PACS numbers: 03.75.Mn, 03.75.Kk 



I. INTRODUCTION 

A Bose-Einstein condensate (BEC) in an optical 
trap exhibits a rich variety of spin-related phenom- 
ena, such as various magnetic phases 2] and spin domain 
formation The nuclear spin of 87 Rb and 23 Na is 3/2, 
and combined with the spin 1/2 of the outermost elec- 
tron, the possible hyperfine spins of these atomic species 
are / = 1 and / = 2. The / = 1 BEC was first realized 
at MIT with 23 Na pj, for which the antiferromagnetic 
behavior was observed Q. On the other hand, 87 Rb in 
the / = 1 hyperfine state was found to be ferromag- 
netic HQ The / = 2 23 Na condensate was also realized 
by the MIT group @, whose ground state phase at zero 
magnetic field is predicted to be antiferromagnetic 0. 
However, the ground state phase and the spin dynamics 
of the / = 2 23 Na BEC have not been studied because 
of its very short lifetime (a few milliseconds), due to the 
fact that the energy of the / = 2 state is higher than 
that of the / = 1 state. The / = 2 87 Rb BEC also lies 
energetically higher than its / = 1 counterpart, but its 
lifetime is much larger (~ 100 ms) due to a fortuitous co- 
incidence of the singlet and triplet scattering lengths . 
The coherent spin dynamics of the f — 2 BEC have been 
observed 0,0. 

The ground state of the / = 2 BEC of 87 Rb is pre- 
dicted to be very close to the phase boundary between 
the antiferromagnetic and cyclic phases 0, EJ • Accord- 
ing to calculations performed by Klausen et al. 0, the 
ground state phase of this BEC at zero magnetic field 
is barely in the antiferromagnetic phase. Recent experi- 
ments performed by the Hamburg group Q and by the 
Gakushuin group [!j appear to support this prediction. 
However, due to the experimental uncertainties, the pos- 
sibility that the ground state phase is cyclic has not been 
excluded. In the Hamburg experiment 0, the m = ±2 
mixture was shown to be stable, which is the ground 
state of the antiferromagnetic phase. However, a mag- 
netic field of 340 mG was applied to the system in this 
experiment. Since the magnetic field lowers the energy 



of the to = ±2 states due to the quadratic Zeeman ef- 
fect, the observed stability of the to = ±2 mixture may 
be due to the presence of the magnetic field. That there 
is no spatial separation between the m = ±2 states is a 
necessary condition for a BEC to be antiferromagnetic, 
but it is not sufficient. In the Gakushuin experiment [jj, 
the spin dynamics starting from the m = state were 
investigated for various values of the magnetic field. It 
was found that the population of the m = ±2 compo- 
nents after 70 ms evolution increased with a decrease of 
the magnetic field, which might imply that the / = 2 
BEC is antiferromagnetic at zero magnetic field. How- 
ever, the initial m — state is a highly excited state and 
the resultant state after 70 ms is far above the ground 
state. Thus, the experimental observations do not ex- 
clude the possibility that the ground state phase at zero 
magnetic field is cyclic. In fact, the Hamburg group ob- 
served that the spin configuration of the cyclic phase has 
a long lifetime (IITH. 

There are two experimental difficulties that hinder the 
determination of the ground state phase. One is the short 
lifetime of < 100 ms of the upper hyperfine manifold of 
87 Rb, which makes the equilibrium spin state hard to 
reach. In order to realize the ground state phase, we 
must equilibrate the system at least for a few seconds, 
as in the / = 1 case 0, ■ The other difficulty is the 
small energy scale (~ 0.1 nK 4J) of the spin-singlet-pair 
term whose sign determines the ground state phase. In 
order to determine the sign of the spin-singlet pair energy, 
the magnetic field must be reduced so that the quadratic 
Zeeman energy is less than the spin-singlet pair energy. 
However, in such a low magnetic field the system may be 
affected by stray ac magnetic fields 0. 

To circumvent these difficulties, we propose a method 
to determine the value of the spin-singlet pair energy by 
spin exchange dynamics. We show that by controlling the 
initial population and phase in each magnetic sublevel, 
we can measure the value of the singlet-pair energy from 
the initial spin dynamics. This method therefore has the 
advantage that equilibrating the system, which takes a 
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long time, is not necessary. Furthermore, the quadratic 
Zeeman energy is allowed to exceed the singlet-pair en- 
ergy and extreme suppression of stray magnetic fields is 
not required. 

This paper is organized as follows. Section [H] formu- 
lates the mean field and Bogoliubov theories of a spin- 2 
BEC in the presence of the quadratic Zeeman effect. Sec- 
tion studies the spin dynamics in a homogeneous sys- 
tem and proposes a method to determine the spin-singlet 
pair energy from the spin exchange dynamics. Section llVl 
numerically verifies the proposed method for a trapped 
system with two-body loss. Section Ivl discusses the ini- 
tial spin preparation and Sec. I VII concludes the paper. 
Detailed derivations of the analytic solutions are given in 
the Appendix. 



II. MEAN FIELD AND BOGOLIUBOV 
ANALYSIS FOR A HOMOGENEOUS SYSTEM 

A. Formulation of the problem 

From rotational symmetry, the low-energy interaction 
between two atoms with hyperfine spin / can be classified 
in terms of the total spin T = 0, 2, • • • , 2/, and each 
scattering channel can be described by the corresponding 
s-wave scattering length a? 0]. I n the case of / = 
2 there are three scattering channels T = 0, 2, 4, and 
the interaction is described by (AttH 2 /M)(aoVo + a 2 V 2 + 
a^P^ )S(ri — r 2 ), where M is the mass of the atom and 
7-y projects the state of two colliding atoms into the 



state of total spin T . The interaction Hamiltonian can 
be rewritten as (c + c\S\ ■ S 2 + c 2 Vo)$(ri — r 2 ), 

where Si and £2 are the spin operators for atoms 1 and 
2 and the interaction coefficients are given by 
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For the linear and quadratic Zeeman effects, the energy 
of the atom depends on the magnetic field B as pm+qm 2 , 
where p = pbB/2 and q — — (/ib-B) 2 / (4Ahf) with p-Q the 
Bohr magneton and Ahf the hyperfine splitting between 
the states / = 1 and / = 2 0|. We assume q < 
throughout this paper, which is the case for spin-2 87 Rb. 

The mean-field energy of the entire system is given by 



where m = — 2, • • • , 2 denotes the magnetic sublevels, V 
is an external potential, 



That 



= E 1^ 



m=-2 

is the total particle-number density, 



is the spin vector with 5x5 spin-2 matrix S, and 
1 
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y/5 



(2^2^-2-2^1^-1+^0) 



(3) 



(4) 



(5) 



is the singlet-pair amplitude. Minimizing the action, 

s= f dt(ihJ2i>*M m -E), (6) 

gives the multicomponent Gross-Pitaevskii (GP) equa- 
tions for the spin-2 BEC as 
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±2 



dt 



= -7T77 V + V" ± 2p + 4g ^±3 + contotV'ia 
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dt 



Cl (F T ^ ±1 ±2F z i; ±2 ) + -±=A 00 r T2 , (7a) 
~2M^ 2 + V±P + q ) ^ J±1 +c ° ntot '^ ±1 



+ci ^^F T Vo + F±V>± 2 ± F,V±i 
c 2 

h 



(7b) 



ih-^r = ( -ttttV^ + y I ^0 + co^totV'o 



dt 



2M 
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+ ^-ci (F_V-i + F+Vi) + ^A 00 r , (7c) 

where F± — F x ± iF y . 

We note that the linear Zeeman term only rotates the 
hyperfine spin about the z axis, and hence the spin dy- 
namics are essentially independent of the linear Zeeman 
term. In fact, setting ij) m — > e lpmt /*"ip m and noting that 
the phase factor also arises in the F± terms as e T * pt / h F±, 
we can completely eliminate the linear Zeeman terms in 
Eqs. H7afl and J7EJ. This property is due to the rotational 
symmetry of our system with respect to the z axis, which 
results in the conservation of the projected angular mo- 
mentum on the z axis, (F z ) — J dr^2 m m\ip m \ 2 . 



E = dr 



S (-^TT^ + V+pm + qmAyj, 



, C 2 



2M 

C-2 



IA,ol 2 



(2) 



B. Ground states in a homogeneous system 

In an ultracold spinor BEC which is isolated from the 
environment, the total spin angular momentum in the di- 
rection of the magnetic field is conserved for a long time 
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(> lsec) @. We therefore minimize the energy of the sys- 
tem in the subspace of a given (F z ) = J dr J2 m m \i ! m^\ 
we shall refer to the resulting minimized state as the 
"ground state" . For simplicity, we restrict ourselves to 
the subspace of (F z ) = in this paper, i.e., we seek the 
ground state in the subspace of (F z ) = 0. 

We consider the uniform case by setting V = and 
ipm = \/nC m , where the density n is constant and Cm 
satisfies ^ m |£m| 2 = 1- The energy per atom is given by 



E 



— = q 2^ m \Cm\ + +y 



C 2, ,2 



(8) 



where / = F/n = J2 m m' Cm S mm'Cm', a 0Q = A 00 /n = 
(2C2C-2 - 2&C-1 + Co 2 )/V5, and 

Cj = Cjn (9) 

for i = 0, 1, 2. The linear Zeeman term is absent in 
Eq. © because (F z ) = 0. 

In order to find the ground state, we compare the en- 
ergies of the stationary states given in Ref. The 
energy of the antiferromagnetic state, 

CAF = (e 2X2 ,0,0,0,e 1 *- 2 ), (10) 

where x-m is the phase of each component, is given by 



£af = 4g + c /2 + c 2 /10. 

This energy is independent of the phase factor e tx " 
the cyclic state, 

Cc = I ~^e iX2 sin 6,0,e tXo cos 0,0, -^=e ix - 2 sin 6 
\V2 V2 

the energy of the system takes the form 



c , C 2 1 2 

ec = y + ^|cos < 



e A sin 



rsm 



(11) 
. For 



(12) 



(13) 



where \ = X2 + X-2 - 2xo- When c 2 < 10|g|, Eq. O 
is minimized by = tt/2, which corresponds to the an- 
tiferromagnetic state £af- When c 2 > 10|g|, Eq. i|T3|) 
becomes minimal with x = 7r and 



2 „ 1 5g 
cos 2 6< = - + — . 
2 c 2 

Substituting this into Eq. I|13|) . we obtain 

10g 2 



(14) 



ec 



2q 



Co 
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(c 2 > 10M). (15) 



The states Caf and Cc satisfy F z (r) = J2 m m\'4>m ( y r )\ 2 = 
0. 

Even if the local spin F z (r) is nonzero, the total spin 
(F z ) can be zero by the formation of the spatial spin 
structure 0|. For example, a staggered domain struc- 
ture of the ferromagnetic state gives (F z ) — 0. The ki- 
netic energy at the domain wall can be neglected when 



the region of the wall is much smaller than the region of 
the domain. Hence, we also consider the case of F z (r) ^ 
0, bearing in mind that the total spin is maintained as 
zero, (F z ) — 0, by some structure formation. The en- 
ergy of the ferromagnetic state, £f = (e* X2 , 0, 0, 0, 0) or 
(0,0,0,0,e iX - 2 ), is given by 

e F = 4g + co/2 + 2c x . (16) 

We also consider a state found in Ref. [TT| . 

Cm = (e iXa cob fl.O.O.e^- 1 sin 9,0) 

or (0,e lxl sin 6,0,0, e lx ~ 2 cos 6). (17) 

The energy of this state is obtained to be 



£ M = q(l + 4cos 2 6) + -j + y (3 cos 2 6 - l) 2 . (18) 

This is minimized by cos 2 6 = 1/3 — q/ (3ci) for ci > \q\/2, 
and the energy then becomes 



£M = 2q+ Y-2c- 1 



(ci > \q\/2). (19) 



Comparing the energies £af, Ef> e C> an d £m ; w e obtain 
the phase diagrams shown in Fig. ^ We have assumed 
here that one of the states Caf, Ccj Cf, and Cm is the 
ground state, based on the numerical results of Ref. [Tlj . 
We have confirmed that this assumption is correct using 
the Monte Carlo method. 

In the case of 87 Rb, ci is positive and|c 2 | is at most 
of the same order of magnitude as c\ [li|. Hence, the 
ground state of the spin-2 87 Rb BEC is antiferromagnetic 
or cyclic depending on whether C2 < 10|g| or £2 > V0\q\. 
Therefore, in order to determine the sign of ci from the 
ground state, we must suppress the magnetic field so that 
the condition \q\ < cWlO is met. In the experiment of 
the Hamburg group 4], a magnetic field of B — 340 
mG was applied, which corresponds to |<?|/fcB — 0.4 nK. 
The mean atomic density was n ~ 4 x 10 14 cm -3 , and 
if \c 2 \M I (4irh 2 ) ~ a B , we have |c 2 |/fc B ~ 1.5 nK. Thus, 
c-2 ^ 1 1 q I , which is consistent with the observed stability 
of the antiferromagnetic state Caf- In the Gakushuin 
experiment B > 100 mG and n ~ 2.1 x 10 14 cm -3 , 
corresponding to |g| > 0.03 nK and \c-Mks — 0.8 nK. 
Hence, the reported parameters of Ref. [|| belong to the 
cyclic phase. 



C. Stability of the stationary states 

The antiferromagnetic (or stretched) state, 



*AF = ^2(1.0,0,0,!), 



(20) 



was found in experiments to be stable in the presence 
of a magnetic field 4]. To examine this result, we will 
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(b) 



20c 




FIG. 1: The ground-state phase diagrams for p — and q < 
(a) with respect to ci and ci and (b) with respect to \q\ and 
c\. The symbols AF, C, F, and M indicate the states £af, 
£c, Cf, and £m, respectively. In (b), Si > is assumed. 



investigate the stability of this state using Bogoliubov 
analysis. 



Substituting Eq. l(2U|l into Eqs. J7J), we find that it 

evolves as e -i^[ e - 2ipt / h , 0, 0, 0, e 2ipt / h ) with fx = Aq + 
co + £2/5, and hence Eq. H20[l is stationary in the rotat- 
ing frame e lmp ip m . We can therefore perform Bogoliubov 
analysis with respect to the state (|20fl in the rotating 
frame, which is described by the GP equations @ with- 
out the linear Zeeman terms. We set the wave function 
as -0 = e~ iflt / h (i& AF + 4>) and substitute it into the GP 
equations in the rotating frame. Taking the linear terms 



0(r, t) = J2 (uke i(k - r -" kt) + vle-^-^A , (22) 

k 

and the eigenenergies are obtained as 



[e k (£ fc + 2c + 2c 2 /5)] 1/2 , 
[e k (e k + 8c 1 -2d 2 /5)} 1/ \ 

3q) (e k -3q + 2d 1 -2d 2 /5)} 1/2 . 



(s k - Aq) (s k -Aq- 25 2 /5)] 



1/2 



(23a) 
(23b) 
(23c) 
(23d) 



where e k = (hk) 2 /(2M). 

If the eigenenergy is complex, the corresponding mode 
is dynamically unstable against exponential growth. The 
eigenvectors Uk and Vf~ of the first mode (|23a|l are pro- 
portional to (1, 0, 0, 0, 1), which is the same form as \I/af, 
indicating that a collapse occurs if 2co + 2c2/5 < 0. 
The eigenvectors of the second mode (|23b|l are propor- 
tional to (1,0,0,0,-1). Since this eigenmode transfers 
the m = ±2 component to the m = =f2 component, 
excitation of the mode (|23b|) gives rise to exchange of 
atoms between the m — ±2 components. This implies 
that phase separation between the two components oc- 
curs for 8ci — 2c2/5 < 0. Therefore, the fact that no 
phase separation is observed in experiments Q only in- 
dicates that C2 < 20ci; it is not sufficient to conclude 
that C2 < 0. The third mode l|23c|l is two- fold degener- 
ate and the eigenvectors are proportional to (0, 1, 0, 0, 0) 
and (0,0,0, 1,0). Therefore, when 2ci - 2c 2 /5 < 3g, the 
state (|20|l is dynamically unstable against the growth of 
the m = ±1 components. Similarly the last mode (|23d|) 
is proportional to (0,0,1,0,0) and describes the growth 
of the m — component. Dynamical instability arises in 
this mode for £2 > 10\q\, consistent with the fact that the 
ground state has the m = component for c 2 > 10|g|, as 
shown by Eq. ((Til . 

For later use, we also perform Bogoliubov analysis for 
the stationary state * = (0,0, -y/n, 0,0). The eigenener- 
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gies are obtained to be 



[e k (e k + 2c + 25 2 /5)} 1/2 , 

[{e k + q) (e k +q + 6c x - 2c 2 /5)] 1/2 

[{e k +4q) (e k +4q-2d 2 /5)] 1/2 . 



(24a) 
(24b) 
(24c) 



The eigenvector of the first mode (|24a|) is proportional 
to (0, 0, 1, 0, 0) oc implying that the state *f? collapses 
if 2c + 2c 2 /5 < 0. We note that this condition for 
the collapse of the BEC is the same as that for *&af- 
The second and third modes, Eqs. l|24bH and H24c|l . are 
both two-fold degenerate. The eigenvectors of the mode 
(|24b)l are proportional to (0,1,0,0,0) and (0,0,0,1,0), 
and therefore the dynamical instability in this mode in- 
creases the m — ±1 components. Similarly, the mode 
(|24c() has eigenvectors proportional to (1,0,0,0,0) and 
(0, 0, 0, 0, 1). In contrast to the excitations for *&af, there 
are always dynamically unstable modes in Eqs. (|24b(l and 
(|24c|) because of q < 0, and therefore the m = state is 
always dynamically unstable for B > in a homogeneous 
system. 



The GP equations then reduce to 

C2 



ih( 



:Co(2C2C-2+Co 2 ), 



(26a) 



C2 



^C±2 = -f C; 2 (2C2C-2 + C ) + 4gC±2. (26b) 



Differentiating Eq. I)26a() with respect to time and using 
Eq. i|26bl) . we can eliminate £± 2 to obtain 



h 2 Co = (e - 4q + 8q\( Q \ 2 ) Co - l + 8q j /,,„. 



where the energy per atom, 



e = ^|2C2C- 2 + Co 2 r + 4g(|C 2 | 2 + IC- 2 | 2 ) 



(27) 



(28) 



is a constant of motion. It is interesting to note that the 
form of Eq. (|27|) coincides with that describing a ID BEC 
on a rotating ring if the time derivative is replaced with 
the spatial derivative 0]. The solution is obtained (see 
Appendix for derivation) as 



III. SPIN DYNAMICS IN A HOMOGENEOUS 
SYSTEM 

A. Analytic solutions 

In order to determine the value of C2 from the spin dy- 
namics, we consider the spin dependent interaction that 
is sensitive only to the spin-singlet term. The elementary 
processes in the spin-singlet channel are + <-> + 0, 
+ 0^1 + (-1), + 0^2 + (-2), and 1 + (-1) <-> 
2 + (—2). On the other hand, the c\ term consists of ele- 
mentary processes such as mi+m 2 <-> (mi + l) + (rn 2 — 1). 
The spin independent interaction does not flip the spin. 
Thus, the elementary process appearing only in the c 2 
term is + <-> 2 + (—2). We therefore focus on this 
process and consider the initial state of the form 



1 



sin6>,0,e IXo cos 0,0, 



1 

7T 



sin( 



(25) 



Since the overall phase and spin rotation about the z 
axis do not affect the physics, the relevant phase appears 
only in the combination of \ 2 + \-2 — 2,Yo- We therefore 
assume that xo = X2 — and < 9 < ir/2 without loss 
of generality. 

If C±i are exactly zero in the initial state, we find from 
Eq. i|7b(l that (±i(t) are always zero within the mean 
field approximation. We therefore assume Chi (t) = 
in the following analysis. In experiments |4|, lll| . this 
assumption holds at least for an initial period of ~ 100 
ms. Although the m = ±1 components exponentially 
grow in the presence of the dynamical instability, the 
assumption C±i{t) = is still valid in the early stage of 
time evolution. 



Co(t) = e^°(*V^osn 2 (af + p \v) + B , (29) 

where <po(t), Ao, Bo, a, (3q, and v are given by Eqs. I|A24|I . 
(|AT6)l . (JXTSJ), l|A2l)l . l|Al2)l . and (|A17|l . respectively. The 
function sn is a Jacobian elliptic function |18| . 
If we assume X-2 = 7r in the initial state, i.e., 



C(o) 



—= sin 6*o, 0, cos 9 , 0, j=sm9 ), (30) 

\/2 y2 



the constants in the solution (|29|l become /?o = 0, Bq = 
Co(0) = cos 2 9q, and 



Aa = ^ 2 (0)[1-Co 2 ( 



(31) 



7 ^{c 2 [l-2C 2 (0)] + 10 (? }, (32) 



5?r 

clCo 2 (0) [l-Co 2 ( 



10q{c 2 [1-2C 2 (0)] + 10 9 }' 



(33) 



The sign of Ao is opposite to that of C2 because q < 
0. The constants a and v are real and positive if 
c 2 [1-2C 2 (0)] + 10g < 0. 



B. Measurement of C2 

We suppose that the initial state 1)30(1 is prepared 
and the magnetic field satisfies \q\ > |c2|/10. Under 
this condition, we find that a and v in Eqs. and 
(|33|l are real and positive. Since the atomic density 
is typically n ~ 10 14 cm~ 3 in experiments 0, Q and 
|(7do — 10a2 + 3a4)/7| is of order of the Bohr radius at 
most 0|, the condition \q\ > | C2 1 /10 requires B > 100 
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20 40 60 20 40 60 

t [ms] t [ms] 

FIG. 2: Time evolution of the relative population of each spin 
component for a density of atoms of 2.1 x 10 14 cm~ 3 9] and a 
magnetic field of 250 mG. The scattering lengths are assumed 
to be (ci, C2)M /(A-Kh 2 ) = (ob, Ob) in (a) and (as, -ob) in (b), 
where as is the Bohr radius. The initial state is chosen to be 
Co = Vol, <2 = -C_2 = v^35, and C±i = \/O001. The 
dashed curves are drawn according to Eq. 1351 . 



also determine the magnitude of ci by measuring the am- 
plitude of the oscillations. In Fig. [21 Eq. ((35(1 is plotted 
as dashed curves, which are in good agreement with the 
numerical solutions (solid curves) in the early stages of 
the time evolution. 

The differences between the numerical and analytic 
curves in Fig. [5] arise mainly from the growth of the 
m = ±1 components in the numerical result due to 
the dynamical instability. We find that the growth of 
the m = ±1 components in Fig. [2] (b) is larger than 
that in Fig. [21 (a). This can qualitatively be under- 
stood from the Bogoliubov energy (|24b|) , which describes 
the growth of the m = ±1 components from the state 
* = (0, 0, y/ri, 0, 0). The imaginary part of Eq. I|24bjl at 
£k = 0, [q(q + 6ci — 2C2/5)] 1 / 2 , is larger for C2 < than 
for C2 > since q < and q + 6ci > in the present case. 

We note that the relative phase x = X2+X-2 — %Xo = t 
in the initial state 1(25(1 is crucial for determining the value 
of C2 by the above method. For example, when the initial 
state is given by Eq. (|2*5|l with x = 0, or without loss of 
generality, 



mG. At this magnetic field, the ground state in the sub- 
space of (F z ) — is the stretched state (120(1 irrespective 
of the sign of C2 . 

The Jacobian elliptic function sa{pd\v) initially in- 
creases for positive a and v, and hence the m — compo- 
nent |Co(*)| 2 m Eq. 1(29(1 initially decreases for C2 > and 
increases for c 2 < 0. Therefore, we can determine the 
sign of C2 from the sign of the initial change in the spin 
population, which is the main result of this paper. This 
method is suitable for atomic species with short lifetimes, 
since only the initial stage of time evolution is needed for 
the determination of C2. 

Figure shows time evolution of the spin populations 
which are obtained by numerically solving Eq. 1(26(1 . The 
results confirm the relation between the sign of C2 and 
the sign of the initial change in the spin populations. In 
Fig. |21 the m = ±1 components are assumed to have 
small initial values (C±i = VO.001) to assess imperfect 
population transfer in realistic situations. We find that 
the dynamical instability in the m = ±1 components 
have little effect on the initial dynamics of the m = 0, 
±2 components. 

In Fig. |21 v ~ 0.01. For < v < 1, the Jacobian el- 
liptic function can be approximated by the trigonometric 
function as 0] 

sn{at\v) = sin erf + 0{v), (34) 
and the spin evolution is then approximated by 

I Co (^) 1 2 - cos 2 O + A Q sin 2 at, (35a) 



IC± 2 (i)| 5 



— (sin 6q — Aq sin at) 



(35b) 



The amplitude of this oscillation is A$ cx — C2, showing 
explicitly that whether the spin population initially in- 
creases or decreases depends on the sign of c-i. We can 



1 

71 



sin 6q, 0, cos 



^o,0, -^=sin6» 



(36) 



the time evolution becomes Eq. ((29(1 . where Ao, v : and 
a are given by Eqs. I(A33() . I(A34(I . and l(A35(l . respec- 
tively (see Appendix for details). These constants satisfi 
C2A0 > and v < 0, and a is real. Using the relation 



lsty 
.18] 



sn(at I — \v\) = 



0" 



sd y/l + \v\at 



1 



and noting that the Jacobian elliptic function sd = sn/dn 
initially increases, we find that | Co (01 2 initially decreases 
for C2 < and increases for ci > 0. This behavior is 
opposite to that in Fig. and therefore control of the 
relative phase x m the initial state is crucial for deter- 
mining the sign of C2. Figure|3shows the initial evolution 
of I Co 1 2 with x varying from to tt, which indicates that 
the phase of the oscillation depends on x- 



C. Dynamics in the cyclic phase 

If |q| < C2/IO, the stretched state 1(20(1 is dynamically 
unstable against growth of the m = component, as 
shown by Eq. ((23d(l . This is because the ground state is 
in the cyclic phase given by Eq. (|T2*)l if \q\ < C2/IO and 
c 2 < 20ci. We can therefore conclude that £2 — 10|q| > 
and hence C2 > 0, if we observe the growth of the m = 
component from the initial stretched state. 

Suppose that cos 2 6 <C 1 in the initial state given by 
Eq. l(3U|) and \q\ < c 2 /10, then v is negative and a is 
pure imaginary from Eqs. I(32|) and l(33|) . In this case, the 
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(a) 0.4 



(b) 0.4 




10 20 

t [ms] 



FIG. 3: Time evolution of |£o| 2 with the same parameters as 
in Figs. [5] (a) and (b). The initial state is given by £q = \/0.3, 
C±i = 0, ^2 = V0.35, and £_2 = exp(j7rj/4)\/0.35, where j is 
an integer ranging from to 4. 



Jacobian elliptic function can be rewritten as 



sn(i|a|f I — |H) 



sd y/l + \v\\a\t 



(38) 

Since \v\ <C 1 and hence 1/(1 + \v\) ~ 1, we can make the 
approximation [18j 



sdHa+nr 1 ) 



sinh it + ^ (sinh u cosh u — u)sechu 



1+4 (sinh tt cosh u — u) tanh u 



(39) 



where it = \a\t. When u ^> 1, Eq. (|39[1 reduces to 

sd(«|(l + M)- 1 ) * -^cch ^|«|* - \ In . (40) 



Equation l|29|) is thus approximated to give 

|Co(i)| 2 =s -^sech 2 C| a |t- iln^ + cos 2 o . (41) 



The dashed curve in Fig. 01 shows Eq. 141J1 . A com- 
parison of this with the corresponding numerical result 
(solid curve) shows that Eq. (|4*T1) is in good agreement 
with the numerical result, with a small discrepancy aris- 
ing mainly from the growth of the m = ±1 compo- 
nents in the numerical result due to the dynamical in- 
stability. The ttl — component reaches a maximum at 




100 200 

f[ms] 



300 



FIG. 4: Time evolution of the spin populations with a mag- 
netic field of 100 mG. The other parameters are the same 
as in Fig. |2] (a) . The initial state is given by fo = V 0.009, 
C±i = VOIM, < 2 = -C-2 = v / 0l>9. The dashed curve is 
drawn according to Eq. (14111 . 



t ~ ln(16/|z/|)/(2|a|) and the maximum value of |Co(^)| 2 
is given by ~ —Aq/\u\ ~ 1 + 10q/c2- 

We note that a characteristic time scale of the dynam- 
ics in Fig. Q]is ~ 100 ms, which is much larger than the 
~ 10 ms seen in Fig. [5] This is because the former time 
scale is set by the dynamical instability, while the latter 
is set by the time evolution from a non-stationary state. 
This method can thus determine the value of C2, provided 
that (i) C2 is positive, (ii) the magnetic field is suppressed 
so that \q\ < C2/IO, and (iii) the time scale of the dynam- 
ics is much faster than the lifetime of the / = 2 state. If 
C2 is negative, the present method cannot determine the 
value of C2 . 

In the case of Fig. 0] (|g| < 62/IO), the m = pop- 
ulation initially increases for C2 > 0, which is opposite 
to the case of Fig. {\q\ > |c 2 |/10). Nevertheless, the 
method to determine the sign of C2 in Sec. 1111 ±31 is still 
applicable, since we can easily find whether we are in a 
region of \q\ > C2/IO or \q\ < C2/IO. If the amplitude 
of the initial oscillation monotonically decreases without 
changing sign with an increase in the magnetic field, we 
can conclude that \q\ > c 2 /10. 



IV. SPIN DYNAMICS IN THE TRAPPED 
SYSTEM 

We investigate the spin dynamics in a trap potential 
to verify the results in Sec. IIIII in a realistic situation. 
We use an axisymmetric trap potential with radial and 
axial frequencies (uj±,oj z ) = 2tt x (237,21) Hz, which 
was used in the experiment of Ref. 0- The two-body 
loss rate was found to be roughly independent of the 
magnetic field and therefore it is considered to be in- 
sensitive to the spin- mixing dynamics 0- We there- 
fore take into account the two-body loss by adding the 
term — ifiK2\'^ m \' 2 "4 ) m/'^ to the right-hand side of Eq. Q 
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FIG. 5: Time evolution of the relative population 
J dr\ip m \ 2 /N of each spin component m in the trap with 
u± — 237 Hz and ui z = 21 Hz [J. The scattering lengths are 
assumed to be (co, ci, C2)M / (A-kTi 2 ) = (100aB, «b, Ob) in (a) 
and (lOOeiB , cib, — Ob) in (b), where ob is the Bohr radius. The 
strength of the magnetic field is 250 mG, and the two-body 
loss coefficient is K 2 = 1.4 x 10~ 13 cm 3 s _1 0. The m = -2 
ground state ?/> 9 is prepared, and is initially transferred to 
each spin component as ipo = V0.298ip g , Tp±i = \/0.001ip g , 



and i>2 



-4>- 



The insets show the column 



density distributions of the m = —2, • • • , 2 components from 
top to bottom at t — 100 ms, each having a 70 x 8/xm field of 
view. 



with K2 = 1.4 x 10 cm s . We prepare the ground 
state ip g for the m = —2 state by the imaginary-time 
propagating method [l9l | and transfer it to each spin 
component as tp = \/0.298ip g , ip±\ = vCTOOT^g, and 

use this state as an initial 
state, where the initial number of atoms is assumed to 
be N(t = 0) = 4.7 x 10 5 @. 

The time evolution of each spin population obtained by 
numerically solving the GP equation JJJ with the two- 
body loss terms is shown in Fig. [S] We find that the 
aforementioned dependence of the initial spin evolution 
on the sign of C2 is valid in the trapped system, i.e., the 
quantity |Co| 2 = / dr\ip \ 2 /N first decreases for c 2 > 
and increases for C2 < 0. The growth in |"0±i| 2 due to 




FIG. 6: Schematic illustration of the energy levels of the 
f — 2 and /' = 2 D\ states. The circularly polarized light 
a± changes m by ±1 and hence starting from the / = 2, 
m — state, the system undergoes transitions only to /' = 2, 
m = ±1 and / = 2, m = ±2 states with the Clebsch-Gordan 
coefficients given by ±1/2 and =Fl/v6, respectively. 



the dynamical instability and the two-body loss does not 
affect this dependence. Thus, we have confirmed that 
the sign of C2 can be determined from the spin dynamics 
in a realistic situation. We find that the growth of the 
m = ±1 components is larger for c 2 < than for c 2 > 0, 
which is similar to the homogeneous case shown in Fig. [21 

The insets in Fig. [S] show the column density distribu- 
tion of each spin component at t = 100 ms. The initial 
distribution is almost preserved at t = 100 ms and no 
spin domains are formed. 



PREPARATION OF THE INITIAL SPIN 
STATE 



As shown in Fig. |3 our method to determine the sign 
of c 2 hinges on the relative phase % 2 + X-2 ~ %Xo i n the 
initial state Therefore we must control the phase 

experimentally. In the experiments of Refs. 0, ||, the 
rapid-adiabatic-passage technique |20j was employed to 
prepare the initial state. This method can control the 
phase in the initial state in principle, if the relative phases 
between the applied rf fields are controlled. We propose 
here a new method to prepare the initial state using the 
Raman transition. 

We consider the Raman transition, for example, be- 
tween the f — 2 and f — 2 D\ states by circularly 
polarized photons, as illustrated in Fig. The initial 
state is assumed to be the m = lower state. Hence, the 
relevant states are the m — 0, ±2 lower states and the 
m = ±1 upper states, whose amplitudes are denoted by 
Coj C±2> an d C±ii respectively. The equations of motion 



9 



are given by 



of each component becomes 



iKC±2 = 
ih( = 



(±2p + 4g)C± 2 T -^e*^'^, (42a) 



(42b) 
(42c) 



where <?± and o>± are the coupling constants and frequen- 
cies of the a± photons, E]j 1 is the transition energy of 
the D\ line, and p' and q' are the linear and quadratic 
Zeeman energies for the upper hyperfine manifold. The 
coupling constants g± are given by the product of the 
amplitude of the circularly polarized light and the dipolc 
matrix element between the s and p states. Using the 
new variables £± 2 = e ±2lpt (±2 and C±i = e8U!± 'C±ii the 
coefficients on the right-hand side of Eq. I|42|) become 
time independent. Assuming ui± = Ed 1 — A ± p + q' 
and A ^> \p — p'\, and eliminating (' ±1 from Eq. 1|42|). we 
obtain 



ih(±2 
ih( 



Aq- 



6A 



C ±2 + ihrt( , 



(43a) 



4A 



(l.9+| 2 + l.9-| 2 )Co-^(^C2-^C-2), 

(43b) 



where iMl = g+g*_/ (2y6 A). If the amplitudes of the a± 
fields are the same, i.e., \g+\ = \g~\, and if the magnetic 
field is applied so that the condition q = — |gt + | 2 /(12A) is 
met, the coefficients of the first terms on the right-hand 
side of Eq. become equal, i.e., Aq — |g T | 2 /(6A) = 

-(I.9+I 2 + l5-| 2 )/( 4A ) = 6 <7- Tne solutions are thus ob- 
tained as 



c 2 = 

Co 

C-2 



_ p -i(2p+6q)t/H _ 



n 



: sin 



-iQqt 



cos v / 2|^j 2 i, 



^-i(-2p+6q)t/h_ 



n* 



VW*, (44a) 
(44b) 

:sm v / 2\n\H. (44c) 



We note that the phase difference is always \ — \2 + 
X-2 — 2xo = ti") which is independent of the phase of 
the circularly polarized light and the applied magnetic 
field. This method is thus suitable for preparing the ini- 
tial state J3DJ. 

We can check whether the prepared state has the cor- 
rect relative phase x = 7r - Rotating the spin state (|25|l by 
7r/2 around the axis on the x-y plane as C — e~ J7rS */ 2 £ 
with Stf, = S x cos 4> + S y sin </>, we find that the population 



IC2I 2 = |C- 2 | 2 = i(3cos 2 



V3 sin 26 cos — cos \ 



- sin 2 6 cos 2 x' 



IC1I 2 

ICol 2 



|C-i| 2 = -sin 2 0sin 2 X ', 

\ / y 

— ^cos 2 9 + V3 sin 26 cos — cos ; 



-3 sin 2 9 cos 2 \ ) , 



(45a) 
(45b) 

(45c) 



where x' = 2<j)+ (x2 — X-2)/2. Let us consider the quan- 
tity 

]l = 2(4|C 2 | 2 + |Ci| 2 ) = l+2cos 2 0-\/3sin26lcos|cosx', 

(46) 

which can be measured by the spin-dependent non- 
destructive imaging reported in Ref. |2l|. We note that 
if X = 7T, Eq. t4lfl) is independent of %' and is con- 
stant, while if x 7^ tt, Eq. (|46|l oscillates, since the phase 
X2 — X-2 rotates in the magnetic field. We can therefore 
confirm that the desired initial state (|30|) is obtained if 
the result of the polarization-dependent phase-contrast 
imaging [2l| is constant in time. 



VI. CONCLUSIONS 

We have proposed a method to experimentally deter- 
mine the sign of the singlet-pair energy of a spin-2 BEC 
from spin exchange dynamics, which determines whether 
the / = 2 87 Rb BEC at zero magnetic field is antiferro- 
magnetic or cyclic. We have obtained analytic solutions 
of the multicomponent GP equations and found that if we 
prepare the initial state in the magnetic sublevels m = 0, 
±2 with the appropriate relative phase relationship, we 
can determine the sign of C2 from the sign of the ini- 
tial change in the spin populations. For example, in the 
case of the initial state in Eq. (|25|l , we can conclude that 
C2 is positive or negative if the m — component first 
decreases or increases. We can also determine the mag- 
nitude of C2 from the amplitude of the oscillation in the 
spin populations. 

Since only the initial evolution for a period of ~ 10 ms 
is needed, we can use this method in the presence of atom 
loss by inelastic collisions and even in the presence of 
dynamical instabilities. We have numerically confirmed 
that this method is applicable for a trapped system with 
a realistic two-body loss coefficient K2- The required 
condition for the magnetic field is \q\ > C2/IO, so our 
method works even when the quadratic Zeeman energy 
exceeds the singlet-pair energy. The initial spin state in 
which the populations and relative phases are controlled 
as Eq. (|25|l can be prepared by, e.g., the Raman technique 
as shown in Sec. Ivl 

We have shown that the coherent spin dynamics serve 
as an efficient probe for investigating spinor properties 
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especially when the atomic species have a short lifetime 
and therefore the equilibrium state cannot be reached. 
The coherent spin dynamics may also reveal other equi- 
librium spinor properties as well as nonequilibrium spinor 
properties that cannot be obtained from the static equi- 
librium. 
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APPENDIX A: DERIVATION OF THE 
SOLUTION (l29l) 



In this appendix, we solve Eq. (12(iH for the initial con- 
dition with xo = Xi — an d < 6 < it/2, and 
derive the solution (|29|l . Defining a new variable 



zo(t) = e^/ 5+4 ^/%(t), 
we can rewrite Eq. I|27|l as 

z = Pz + Q\z \ 2 z , 

where 



1 /25 2 



,~2 



P = 77 "7^ ~^-l&q 



16c 2 q 



25 



Q 



16c 2 q 

5h 2 ' 



(Al) 
(A2) 

(A3) 
(A4) 



with the spin-dependent energy 



£s = ^|2C 2 C- 2 + Co 2 | 2 +4<z(l-|Co| 2 ). (A5) 

The initial conditions for zq are zo(0) = Co(0) = cos#n 
and 



io(0) 



h \ 5 

ic 2 



Co(0) 



^Co(0)[2e^-C 2 (0) 2 +Co 2 (0)], (A6) 



where we have used Eq. (|26a|) . 

We write the complex variable as zo(t) = R(t)e l< ^^ 
with real functions R(t) and 4>(t) : where (f)(0) = and 
z (0) = R(0). Substituting this into Eq. (|A2I) and taking 
the real and imaginary parts, we obtain 



R- PR- QR - Rcj) 



;,2 



0. 



2R<P + R<P = 0. 



(A7) 
(A8) 



These equations can be integrated to give 

r 2 
Z± 
R 2 ' 



R 2 



PR 2 + |i? 4 + C R 

^2' 



(A9) 
(A10) 



where Cr and dp are constants of integration. The time 
derivative of zq at t = 0, i (0) = R(0) + i<fi(0)R(0), gives 
Rei o (0) = R(0) and lmi (0) = <p(0)R(0). The constant 
C,), is then written as 

C+ = 0(O)i? 2 (O) = z (0)lmi (0). (All) 
The constant Cr is similarly obtained as 

^ = |z (0)| 2 -Pz 2 (0)-|z 4 (0). (A12) 

We introduce a new variable p through 

R 2 = A aP 2 + B Ql (A13) 

where Aq and Bq are constants. Substituting Eq. (|A13|) 
into Eq. I|A9(1 . we obtain 

3Q 
2 



^-B 2 + 2PB t ,-< r C l; ) .1,, 
Q t->3 i dd2 i /-»_ d_ \ „— 2 



-Bq + PBZ + C R B -Ci\p 



(A14) 



We determine Bq so that the p 2 term in Eq. (| A 14(1 van- 
ishes: 

^Bl + PBq + CrBq - C\ = 0. (A15) 

The right-hand side (rhs) of Eq. ilA14(l then becomes 
quartic with respect to p. We determine Aq so that the 
rhsofEq. (| A 14(1 becomes proportional to (1— p 2 )(l — vp 2 )\ 
the result is 



An 



1 r 
2Q 



-3QBq - 2P 



±^-(3QB Q - 2P){QBq + 2P) - 8QC R 



with 



QA 



QA + 2P + 3QBq 
Equation (|A14I) thus becomes 
.9 Q A o, 



2v 



-{\-p 2 )(\-vp 2 ). 



(A16) 



(A17) 



(A18) 
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This equation can be integrated to give 
rp(t) 



dp 



We rewrite this equation as 
dp 



where 



V(l-P 2 )(l~^ 2 ) 



0o = 



QA 
2v ' 
p(o) 



2v 



at + (3 , 



dp 



(A19) 
(A20) 

(A21) 
(A22) 



Equation (|A26(1 has the same form as Eq. (| A2|> , in which 
P and Q are replaced by P' and Q', and therefore the 
derivation from Eq. I|A7(I to Eq. (| A24|l can similarly be 
applied. 

We consider the case of X-i = 71 m the initial condition 
given in Eq. I|25|) . In this case, Eq. {2B 

is pure imaginary, 

and Eq. (|A15|) c an be solved to give B = Zq(0) = Co(0)- 
Equation (|A16|I then reduces to 



An 



1-2C 2 (0) 



Wq 

C2 



C2 

Wq 



C 2 (0) [i-Co 2 ( 



(A28) 



where the upper and lower expressions correspond to the 
plus and minus signs of ± in Eq. (|A16|) . Using this Ao, 
v and a become 



The left-hand side of Eq. (|A20I) is an elliptic integral of 
the first kind and hence p(t) — sn(at + PoW) [13 • Thus, 



I Co (*) I = R(t) = VA srf(at + f3 Q \v)+B . (A23) 

From Eqs. IjAlfl and IjAlOfl . the phase cpa of Co is given 
by 



1 f~C2 



dr- 



C,:, 



l A Q sn 2 (at + o \v) + B ' 
(A24) 

Similarly, we can obtain the solution for C±2(t). From 
the form of Eq. I|2tj[l and the initial condition C-2(0) = 



e' x - 2 C2(0), we find that the relation C- 2 (i) = e lx - 2 C, 2 {t) 
always holds. tlence, we consider only £ 2 (i). Eliminating 
Co(t) from Eq. I|26|l . we have the equation of motion for 
C 2 (i) as 



10g{c 2 [1-2C 2 (0)] +10g} 

e;Co 2 (o) [i-Co 2 (o) 
siC 2 (o) [i-Co 2 (o) 



Wq{c 2 [1-2C 2 (0)] +10g}' 

S^ Co 2( 0) [i-c 2 (0)], 



^{c 2 [1-2C 2 (0)] + 10 9 }. 



(A29) 



(A30) 



It follows from Eq. jXT3j> and B = z 2 (0) that p(0) = 
and hence (3q = 0. Using the relation 



sn(ailf) = —=sn{^/vat\v ), 
\ v 



(A31) 



h 2 ( 2 



^(e s +4q)-16q 2 
5 



32c 2 q l/ . |2 



-IC2 



2c 2 • 

C2 -l — K2- 



The new variable z 2 {t) 
to 



(A25) 

js 2 t/(5h)^ reduces Eq . JX251) 



z 2 =P'z 2 + Q'\z 2 \ 2 z 2 , 



(A26) 



where P' = (2£ 2 e s /5 - c|/25 - 16q 2 + 8d 2 q/5)/h 2 and 
Q' = —32c 2 q/(5h 2 ). The initial conditions for z 2 are 
22(0) =Ca(0) = sm6 /V2 and 



i 2 (0) = ^C 2 (0) [1-2C 2 (0) 



5h 



Aiq 



we can show that the two expressions in Eqs. JA28J), 
(|A29|I . and l|A30|l are equivalent to each other. We used 
the lower expressions in Sec. IIIII The m = ±2 compo- 
nents C, 2 {t) = — C- 2 (t) have the same form as Eq. I|A23(I . 
in which Aq, (3q, and Bq are replaced by A 2 , /3 2 , and Bq. 
These constants are given by 



A 2 = 



2C 2 2 (0) 



bq 

c 2 



-^Cl(o) [1-2C 2 (0)], 



C 2 (o). 

(A27) 02 = 0, and S 2 =C 2 2 (0)- 



(A32) 



Finally, we consider the case of X—2 = in the initial condition (|25|l . i.e., all the components are real and positive. 
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We can solve Eq. I|A15(1 also in this case to give Bq — z 2 (0) and hence (3q — 0. The constants Aq, v, and a reduce to 

, {c 2 [l - 2C 2 (0)] + 10q} - \q\y/(c2 + 10 9 ) 2 -40c 2 gC 2 (0) 



2c 2 g I. -V j 

(A33) 

g {g 2 [l - 2C 2 (0)] + 10g} - Mvg + ^g) 2 ~ 40c 2 gCg(0) 34 
9 {C2[1 - 2C 2 (0)] + 10 9 } + |g|V(52 + 10g) 2 - 40^(0) ' 

- - JW- < A35 » 

where we take the lower sign in Eq. (|A"TB|) . Since {c 2 [l - 2Co(0)] + I0q} 2 < (c 2 + 10q) 2 - 40c 2 gCo(0), we nnd Cz A > 
and ^ < 0. , 
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